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Abstract 

The Walker et al. (1989) self-consistent method to predict both the 
elastic and inelastic effective material properties of composites is 
examined and compared with the results of other self-consistent and 
elasticity based solutions. The elastic part of their method is shown to 
be identical to other self-consistent methods for non-dilute reinforced 
composite materials; they sire the Hill (1965), Budiansky (1965) and 
Nemat-Nasser et al. (1982) derivations. A simplified form of the 
non-dilute self-consistent method is also derived. 

The predicted, effective, elastic, material properties for a fiber 
reinforced material using the Walker method was found to deviate from the 
elasticity solution for the U 31 , K 12 . and P 31 material properties (fiber 

is in the 3 direction) especially at the larger volume fractions. Also, 
the prediction for the transverse shear modulus, Pi 2 > exceeds one of the 
accepted Hashin bounds. Only the longitudinal elastic modulus E 33 agrees 
with the elasticity solution. The differences between the Walker and 
elasticity solutions are primarily due to the assumption used in the 
derivation of the self-consistent method, i.e. the strain fields in the 
inclusions and the matrix are assumed to remain constant, which is not a 
correct assumption for a high concentration of inclusions. 


♦Summer Faculty Fellow at NASA Lewis Research Center. 


1 



A summary of some other methods which predict the effective material 
properties of a composite is also presented. In particular, the two 
different bounding solutions of Hashin (1965, 1972) for the transverse 

shear modulus, p 12 , are reviewed. Other methods examined include the 
dilute self-consistent method by Eshelby (1957) and Russel and Acrivos 
(1972), the Mori-Tanaka method (1973), and the three-phase model of 
Christensen and Lo (1979b). 


Introduction 

Over the last 30 years there has been a tremendous effort to predict 
the effective material properties of a composite material based on the 
material properties of the individual constituents. The motivation behind 
this effort is to better design and analyze two part composite materials 
without actually having to produce and test each and every possible 
composite material. Recently, Walker et al. (1989) have proposed using a 
self-consistent method to predict the elastic and inelastic material 
properties and response. This method was developed for Implementation in a 
finite element code to predict the response of a composite structure. This 
report will examine in detail the elastic portion of this solution, and 
will compare it to other self-consistent methods and elasticity solutions 
to assess its ability to predict the effective material properties of a 
reinforced composite material. 

Some of the original work in developing the self-consistent concept 
was done by Hershey (1954), Kroner (1958), and Kerner (1956) on single 
crystal and polycrystalline materials. This work, however, did not examine 
multiple phase materials such as spherical or cylindrical fiber reinforced 
composite materials. In 1957, Eshelby introduced his transformation 
tensor, known as the Eshelby Tensor, which relates the strain of an 
inclusion constrained inside an Infinite elastic matrix to the strain of 
the same inclusion when placed outside of the matrix without any 
constraints. This basic tensor has far reaching implications, one of which 
allows the calculation of the effective material properties of a two phase 
material. One main assumption in the derivation of the Eshelby Tensor is 
that multi-inclusions do not effect one another. For low concentrations, 
this method - known as the self-consistent method for dilute concentration 
(DSC) - accurately predicts the effective material properties. 

Due to the limitation of the DSC method to only low concentrations of 
inclusions, Hill (1965) and Budiansky (1965) simultaneously developed a 
self-consistent method that extends the basic Eshelby method to high 
concentrations - so called the self-consistent method for non-dilute 
concentration (NSC). This method basically lets the matrix material assume 
the properties of the desired effective material properties (still 
unknown), and then places a single inclusion Into this matrix. The Eshelby 
tensor then figures into the description of the strain in both the 
inclusions and the overall composite. The effect of multiple inclusions is 
assumed to be modeled by allowing the matrix to have the effective 
material properties. However, the strain fields in both the Inclusion and 
matrix are still assumed to be constant, as required when using the 
Eshelby Tensor, which is not valid for densely packed inclusions. The 
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prediction of elastic properties using NSC method for material represen- 
tation is better then that of the DSC method. Predictions for the NSC 
method approach the limiting value of the inclusion properties as the 
inclusion concentration ratio approaches one; predictions for the DSC 
method do not. It should also be noted that the solution of the NSC 
problem is iterative, since the effective material properties are imbedded 
directly In the formulation. 

This paper Is mainly concerned with the NSC method, since the 
solution of Walker et al. (1989) is sin extension of this method into the 
Inelastic range. The basic elastic portion of this method is shown to be 
identical in form to the NSC method. Furthermore, the Nemat-Nasser and 
Taya (1981) derivation of the self-consistent method is also shown to be 
the same as the NSC method. 

Another method that uses the Eshelby Tensor is the Mori-Tanaka (1973) 
method. This method accounts for the Interaction of the inclusions by 
defining a second transformation matrix (or strain concentration tensor) 
between the actual matrix material and the inclusions (recall that the 
original transformation matrix describes the relation between an inclusion 
and the effective material properties or homogeneous composite material). 
By including this additional constraint, the effective material properties 
will match the limiting values at both low and high concentrations; plus, 
all five of the elastic constants agree with the elasticity solutions for 
fibrous composites. However, this method is currently limited to elastic 
solutions only. 

In addition to the self-consistent methods for predicting the 
effective material properties, there has been a number of elasticity based 
solutions. Hashin and Shtrikman (1962a, 1962b, 1963) laid much of the 
basic foundation for employing the elasticity equations using energy 
methods to develop solutions to the two part spherical and cylindrical 
inclusion problems. Later, Hashin (1965, 1972, 1979), Hashin and Rosen 
(1964), and Whitney and Riley (1966) presented exact solutions for most of 
the material properties. The one exception is the shear modulus of a 
composite with spherical inclusions, or the transverse shear modulus of a 
fiber reinforced composite. In these cases, only bounds are obtained. The 
Hashin bounds are used extensively in the literature to determine what is 
acceptable, but there is considerable confusion on how to get these bounds 
and what bounds to use. These bounds are discussed In detail later. 

This paper briefly presents each of the basic methods mentioned above 
so that they can be used to compare results with those of the Walker 
method. The final section presents Walker’s method, and derives many of 
the basic equations to show its similarity to the other self-consistent 
methods. A comparison of the various methods is included. 

A quick note on the notation. Bold letters represent tensor 
quantities. The subscripts m and f represent the matrix and fiber (either 
cylindrical, ellipsoidal, or spherical particulate), respectively. This 
report assumes the fibers are always oriented in the 3 direction, such 
that E 33 is the modulus in the fiber direction. The other two directions, 
1 and 2, are perpendicular to the fiber and to each other. Since only 
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transversely Isotropic materials are considered for fiber reinforcing - 
i.e. the fiber array Is hexagonal - the 1 and 2 directions are inter- 
changeable. The five Independent material properties used In this report 
for fiber reinforced composites are: the longitudinal Young modulus, E 33 , 
the longitudinal Poisson ratio, v 31 , the plain strain transverse bulk 
modulus, Ki 2 » the longitudinal shear modulus, P 31 , and the transverse 
shear modulus, p 12 . 


Self-Consistent Method for Dilute Concentrations (DSC) 

All self-consistent methods rely on two conclusions reached by 
Eshelby (1957). First, the free strain, c T , of an Inclusion is linearly 
related by the Eshelby Tensor to the constrained strain, e c , of the same 
inclusion imbedded in an elastic body. Second, the strain is uniform in 
the inclusion. These two conditions allow the effective material 
properties for dilute concentrations to be determined. 

Since the Eshelby tensor is central to the DSC method, it will now be 
explained using the formulation of Russel (1972, 1973), Christensen 

(1979a) and Benveniste (1987). Assume that an ellipsoidal shaped body 
(later to be defined as the inclusion) undergoes a free strain, e T , 
without any constraints, i.e. <r* j = 0. Next, force this inclusion body 
back to its original ellipsoidal shape by applying tractions at its 
surface. Now place this ellipsoidal body inside an Infinite elastic body 
that has an elliptical shaped void of the same shape as the original shape 
of the inclusion. When the inclusion is first placed in the void, just 
before the tractions are released, the infinite body is strain and stress 
free. Finally, the tractions on the Inclusion are released, causing both 
the inclusion and body to change shape. At the interface, both the 
inclusion and the body will have the constraint strain state, e c . The 
Eshelby tensor simply relates the free strain, € T , to this constraint 
strain, e c , as 

e c = Sg t , (1) 

where S is the Eshelby tensor. The S tensor Is a function of the Inclusion 
shape (must be ellipsoidal In shape) and the material properties of the 
body. Exact solutions for S have been worked out for various inclusion 
shapes for isotropic and transversely Isotropic materials by Eshelby 
(1957), Mura (1987), Christensen (1979a), Walker (1989), and others. 

The Eshelby Tensor can be used In composite analyses where the 
inclusion has material properties different from those of the constraining 
body, and where the complete composite undergoes an overall or applied 
strain, g°. This is done by modeling the free strain effect as a material 
property effect. The total strain, e, for an inclusion of the same 
material as the matrix in which It is embedded equals the overall strain, 
e°, plus the constraint strain, g c , less the original free strain, e T , 
giving 


D ffl g (=<r mi ) = D m (c° + c c - g t ), 


( 2 ) 


4 


where D m is the material stiffness matrix for the matrix material. For an 
inclusion of a different material, say D f , the strain does not have an 
original free strain associated with it, since it is modeled as a material 
property effect. Therefore, the total stress is given by 

Df e (so* 1 ) = D f (c° + e c ) . (3) 

Since the stress states for both matrix type Inclusion, c 4 ” 1 , and for the 
fiber type inclusion, must be equal for the same inclusion geometry, 

Eqs. 2 and 3 can be combined to give 

D m ( e o + C c _ e T) = Df (e° + e c ) , (4) 

or by using the definition of the inclusion strain, e f = e° + e c , Eq. 4 
can be rewritten as 

e T = Djj; 1 (D m - Df) e<\ (5) 

Substituting Eq. 5 into Eq. 1 gives 

e c = S m [D^ 1 (D m - D r )] e<\ (6) 

where the components of the Eshelby tensor, S m , are dependent upon the 
matrix material’s properties. Using the definition of the Inclusion strain 
once again gives 


e f = e° + S m [D ; 1 (D m - Df)) e f . (V) 
Rearranging Eq. 7 gives 

e f = Tc° = [I + S m [D * 1 (Df - DJ ]]" 1 e°, (8) 

where T is called the strain concentration tensor, and I is the identity 
matrix. A second equation is now needed. One that relates the effective 
material properties, D, to the material properties of the matrix and 
fiber, and to the strain concentration tensor, T. 

Following the derivation of Russel and Acrivos (1972), the volume 
averaged stress in the inclusions and matrix is added together to give the 
volume averaged stress in an arbitrary subvolume, i.e. 


«»-u VT 

1 N 

0 4 * dV + ry V 
h V L 

| <r f dv. 

(9) 

V-£ 

V n n=i 

= 1 

J v n 



where V is the total volume, V n is the volume of each inclusion, and V-£V n 
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is the matrix material volume. Also, <r® and <r f represent the stress in the 
matrix and inclusion (or fiber) material, respectively. The first integral 
is split into two parts, one part representing the total subvolume as a 
matrix only material, i.e. no inclusions, and a second part representing 
the inclusions as matrix material. The stresses in each of the three 
integral terms Eire then written in terms of strain, 

<cr> = ~ f D m c° dv - ^ J f D m (c° + e c ) dv + 

J V n= i V n 

1 N 

+ y £ J D f (c° + e c ) dv. (10) 

n =i J V n 

The strain in the first integral represents the actual strain of the 
matrix material for the total subvolume as if there are no inclusions. The 
inclusions are modeled by the second and third integrals. The second 
integral subtracts out the inclusion volume made of matrix material, D m , 
whereas the third integral adds back in the actual inclusion material, Df. 
In both cases the strain includes the actual strain, e°, and the 
constraint strain, e c that any Inclusion will experience. Equation 10 can 
be further simplified to 

1 N 

«r> = D m e° + i £ j (D f - D m ) dv. (11) 

n=l J V n 

The effective material property matrix, D, relates the average stress and 
strain over a representative subvolume, i.e. <<r> = D <e>. Assuming that 
the inclusion strain e f does not vary over the volume of the inclusion, an 
assumption that is valid for a single ellipsoidal inclusion, Eshelby 
(1957), Eq. 11 can be integrated to give 

D <e> = D m c° + Vf (Df - D m )e f , (12) 

where V f <= [0,1] is the volume fraction of the Inclusions. Noting that <e> 
and e° are the same far from the inclusion, Eq. 12 can be rearranged to 
give 


D = D m + Vf (Df - D J T, (13) 

where e f = Te°. 


Combining Eqs. 8 and 13 gives the final form for the effective elastic 
moduli of the DSC method, I.e. 

D = D m + Vf (Df - DJ [I + S m [D^ 1 (D f - DJ)]" 1 . (14) 
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This equation can be used for both spherical or cylindrical Inclusions, 
the only difference will be In the values for the components of the 
Eshelby tensor, S m . 

The five effective material properties for a transversely isotropic 
fiber reinforced composite predicted by Eq. 14 are plotted in Figs. 1-5. 
In four of the five material properties, the DSC deviates from the 
elasticity solutions or bounds because of the assumption of no Inter- 
actions between inclusions. For dilute concentrations, the agreement is 
reasonable, and this method can be used. The shear modulus and Poisson 
ratio used for the glass fiber and polyester matrix material properties in 
the Figs. 1-5 are, G m = 0.593 GPa and G f = 29.1 GPa, and = 0.45 and 
v c = 0.21 (Richard (1975)). 

The Eshelby tensor for a transversely isotropic, continuous, 
cylindrical fiber reinforced, composite material - such as the one 
described in the figures - is given by (Walker et al. (1989)) 
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where u m is the Poisson ratio of the matrix. It should be noted that the 
three terms in the lower right quadrant are a factor of two larger then 
the values reported by Eshelby, 1957 and others due to notation. If Voigt 
notation is used, i.e. Sjj, the above values are correct; they are a 
factor of two larger then those associated with the 4 th order tensor 
notation, i.e. S IJkl . A different Eshelby tensor is required for a 
spherical, reinforced, composite material. 


Mori -Tanaka Method 

The Mori-Tanaka method is similar to the DSC method, except it 
enforces the condition that the effective material properties must match 
the properties of the inclusion material as the concentration approaches 
one, i.e. D -> D f as V f -» 1. This is accomplished by introducing a new 
strain concentration tensor, G, which relates the strain in the matrix to 
the strain in the inclusion, i.e. 

e f = Ge ro , 


7 


and by rewriting Eq. 13 as 


D = Dm + Vf (Df - D m ) A, (16) 

where e f = Ae°. 


The A matrix is now assumed to include the effects of multiple inclusion 
interactions, as required for the non-dilute case. Recall, the T matrix 
only assumed a single Inclusion in a matrix with no considerations of 
multiple inclusion interaction. Combining the definitions of A and G, and 
eliminating e f , gives 

Ge m = Ac 0 . (17) 


Because the overall strain, c°, is simply the sum of the relative volume 
fractions of strains in the matrix and the inclusions, the equation 
becomes 


Ge m = A ( V m c m + Vf€ f ) , 


(18) 


where V ro and V f are the matrix and fiber volume ratios, respectively, such 
that V ro + Vf = 1. Rearranging Eq. 18, and using Eq. 15 gives 

A = G [ V m I + VfG]" 1 . (19) 

Equation 19 defines the new strain concentration tensor between the 
fiber strain e f and the applied strain e°. However, G is still not known, 
Mori and Tanaka (1973) make the assumption that at high concentrations the 
matrix # strain equals the applied strain which lead them to postulate that 
G = T, Making this substitution into Eq. 19, and then substituting back 
into Eq. 16 gives 

D = D m + Vf (Df - DJ T [ V ro I + VfT]" 1 , (20) 

where T = [i + S m [D^ 1 (D f - DJ]]" 1 . 

This is the final form of the Mori-Tanaka method. It is similar in form to 
the DSC method, except for the additional term [V m I + V f T] _1 which acts as 
a correction factor at non-dilute conditions. At the low and high limits 
of inclusion concentration, Eq. 20 becomes, respectively, 

D D m as Vf -> 0, and D -> D f as V f -> 1. 

Notice that there is an inconsistency in this derivation; the assumption 
that e m = c°, implying that G = T, when going from Eq. 19 to 20 cannot be 


The theoretical basis for this assumption is discussed in further detail 
in the "Revised Mori-Tanaka Using the NSC Method " section of this report. 


- 8 - 



applied earlier in the derivation, say at Eq. 17. Additional comments on 
the Mor 1 -Tanaka method can also be found in the review papers by 
Benveniste (1987) and Christensen (1990). 

In a more recent paper, Benveniste (1990) has shown that the methods 
of Willis (1981) and Levin (1976) are equivalent to that of Mori-Tanaka 
described above. It is worth noting, however, that the P matrix used by 
Benveniste, defined as P ■ is not symmetric as he states. However, 
the final effective stiffness matrix D is symmetric, as required, even if 
P is not symmetric. 


The five material properties predictions for a transversely isotropic 
composite using the Mori-Tanaka method, Eq. 20, are plotted in Figs. 1-5 
along with other methods that will be discussed in the following sections. 
The Mori-Tanaka results match the elasticity solutions for 4 of the 5 
material properties (Figs 2-5), while the fifth property, the transverse 
shear modulus, matches the lower Hashin Bound (Fig. 1). The material 
properties and the Eshelby tensor are the same as those which were used in 
the DSC predictions. One distinct advantage of this method is that the 
effective material properties can be determined without an iteration 
process, which is not true for all methods as will be seen in later 
sections. Furthermore, the results are in agreement with both the known 
elasticity solutions and the Hashin bounds, but yet can be written in 
compact tensor equation form. 


Hashin Bounds 

One of the first efforts to predict the effective material properties 
of a two phase composite material, for any concentration, using varia 
tional energy principles was done by Hashin (1962), which was later 
generalized to a multi-phase material by Hashin and Shtrikman ( 1962a, b and 
1963). Note that the Hashin method is not a self-consistent method, but an 
exact elasticity method. By applying the variational theorems of the 
theory of elasticity, they determined upper and lower bounds on the 
effective bulk and shear properties for composites with spherical 
inclusions. Basically, first they used a general elasticity formulation to 
obtain the displacement fields of spherical bodies. Then they applied two 
different boundary conditions; one for strain, and the other for stress. 
Finally, the theorem of minimum complementary and potential energies were 
applied to the strain and stress boundary conditions, respectively. Since 
two different boundary conditions were used, two bounds were determined, 
an upper bound and a lower bound. For a spherically reinforced composite 
material, they found the bulk modulus bounds coincided, giving an exact 
solution, but the bounds for the shear modulus did not coincide. 

These same methods were later applied to fiber reinforced, trans- 
versely isotropic, composite materials by Hashin (1964). He was able to 
show that for four of the five properties (i.e. the longitudinal elastic 
modulus, E 33 , the longitudinal Poisson ratio, i> 3 i, the plain strain bulk 
modulus, K 12 , and the longitudinal shear modulus, M 31 K both the upper and 
lower bounds, due to the minimum potential and complementary energies, 
coincide giving an exact solution. However just like the sphere model, the 
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upper and lower bounds for the transverse shear modulus, ji 12 , do not 
match. Because there Is not an exact solution for /i 12 , a number of other 
methods have been developed over the years Including the Self-Consistent 
(Hill, 1965 and Budiansky, 1965), Generalized Self-Consistent (Christensen 
and Lo, 1979), Mori-Tanaka (Mori and Tanaka, 1973), Modified Mori-Tanaka 
(Luo and Weng, 1987), and Differential methods (Roscoe 1952). Although 
these other methods can also predict all material properties, this report 
concentrates on the ambiguous p 12 term. 

Because of the strong theoretical basis of the Hashin bounds, all 
other methods are generally compared with it as a check for correctness. 
If a particular method does not match the Hashin bounds for the four known 
properties, or falls outside of the upper or lower bounds for the fifth 
property, p 12 , then the method is considered to be in error. However upon 
reviewing the literature, one can very quickly become confused with the 
many different bounds that sire referred to as the ’Hashin Bounds’. The 
main cause for this confusion Is that Hashin actually derived two sets of 
bounds for a fiber reinforced composite with randomly placed parallel 
fibers of different sizes. Hashin (1979) also developed a third set of 
bounds for a composite constructed with identical fibers arranged in a 
hexagonal pattern which Is not discussed in this paper. 

The original Hashin bounds for p 12 that were set forth In his 1965 
paper required the solution of 8 linear equations. However, if the fibers 
are assumed to be solid (not hollow) the number of equations reduce to 6 , 
as was presented by Hashin (1972) in a later paper. The upper ( + ) and 
lower (-) bounds, based on the potential and complementary energies, 
respectively, are 


M12 - Mm 
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where v m = Poisson’s ratio of matrix, 

= shear modulus of matrix, 

V f = fiber volume fraction, and 

A 4 = constant from the solution of the 6 linear equations. 


The 6 linear equations, in matrix form, that must be solved to determine 
the constants A 4 are 
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These two bounds, referred to as the composite cylinder assemblage 
(CCA) bounds by Hashin, are plotted in Fig. 6 as a function of the fiber 
volume fraction. In later papers, Hashin (1972, 1979) presented the 
following closed form solutions to Eqs. 21 and 22 for the 6 linear 
equations; they are 
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Kn, and K f = transverse bulk moduli for the matrix and fiber. 


There has been considerable confusion in the literature regarding the 
CCA bounds, because of the two different forms of these equations. Unless 
one reads Hashin’s papers carefully, it is easy to assume that Eqs. 21 and 
22, and Eqs. 25 and 26 are two different sets of bounds; they are not. 
Furthermore, if one uses the standard definition for the bulk modulus, 

K = 5^1 -2 l>) or K = X + 2 P /3 * (27 ) 


for either the matrix or fiber material, Eqs. 21 and 22 will give 
different results from those of Eqs. 25 and 26. Hashin, as well as others, 
e.g. Christensen (1979a), define K* as the plain strain or transverse bulk 
modulus, which is different than the standard isotropic bulk modulus, K. 
The K* term originated from the transversely isotropic, engineering, 
material property K 12 , which is defined as K 1Z = £(C n + C 12 ) where C tJ 
are the stiffness moduli. This transversely isotropic definition is then 
carried over to describe the bulk modulus of both the isotropic fiber and 
matrix materials. Therefore, the correct definition to use with Eqs. 25 
and 26 is 
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(28) 


To make matters more confusing, Hashin (1965) developed a second set 
of upper and lower bounds for ji 12 using the Polarization Extremum 
Principle that was developed by Hashin and Shtrikman ( 1962a, b) and Hill 
(1965) for spherical inclusions. A detailed derivation can also be found 
in Hashin (1972). The final upper and lower bounds for the transverse 
shear modulus of this method are 
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where: 



Vf 


jif = shear moduli of the matrix and fiber, 

Kf = transverse bulk moduli of the matrix and fiber, and 
= fiber volume fraction. 


These equations, referred to as the Arbitrary Cylindrical Phase (ACP) 
model, are plotted in Fig. 6 along with the non-dilute self-consistent 
method (NSC). Note that they do not coincide with the previous bounds, 
Eqs. 21 and 22, or Eqs. 25 and 26. Because the two sets of bounds - the 
CCA and ACP - do not coincide, they tend to be confused with one another. 
In fact, both bounds are still being used and both are labeled as the 
’Hashin Bounds’. Hashin (1979) summarized both sets of bounds briefly and 
stated that the upper bound of the ACP method along with the lower bound 
of the CCA method should be used. This gives the tightest bounds for the 
unknown transverse shear modulus Pi2> and are the Hashin bounds plotted in 
Fig. 1. 


Three-Phase Model 

The Hashin model discussed above gave the exact solution for 4 of the 
5 effective material properties, but only the bounds for the remaining 
property, the transverse shear modulus. Christensen and Lo (1979) later 
developed an exact solution for the transverse shear modulus by using a 
three phase geometry. Like the Hashin solution, they use the basic 
elasticity solution for cylindrical geometries, and thus eliminate the 
need to use the Eshelby solution for Inclusions. Although the three-phase 
model is not a self-consistent model, it is often referred to as the 
generalized self-consistent model because it uses three phases, including 
an outer phase that is modeled as an equivalent, homogeneous or smeared, 
composite material similar to the other self-consistent methods. 

The geometry used in this model consists of a cylindrical fiber 
embedded in a larger cylindrical matrix, which is placed in an infinite, 
homogeneous, composite material that represents the effective composite 
material. The method starts by assuming the solution for the displacements 
in all three phases, and then the eight boundary conditions are applied. 
The displacements, u r and u e , in the outer phase are 
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U, 


_b 

r « " 4m 


2r . , _,b A fb 

— + a 3 (r )+!)- + c 3 |p 


cos (20) , 


u e * = 


4 li 


-2 r , . ,b fb 

— - a 3(»-D F * =3 p 


sln(20) , 


where b is the radius to the outer edge of the matrix material 
In the matrix phase, the displacements are 


Urm " 4p„ 


a 2 (7) ra -3)^ + d 2 £ + c 2 (T? m +l)p + b 2 ^j 


cos (20) 


uem = 


4p m 


r i 3 

a 2 (T) m -3) £ - d 2 P 


sin(20) 


where a s r s b. And in the fiber phase, they are given by 


u rf = 


4/Jf 


ai(r)f-3)|g + d t £ 


cos (20) , 


U rf = 


4Pf 


a!(T}f+3) 


GO 


dl H 


sin(20) , 


where r s a. The constants a 3 , c 3 , a 2 , d 2 , c 2 , b 2 , a*, and d 
determined, and 7) = 3-4y, T) m = 3-4i> m , and T)f = 3-4i> f - 

In addition to the eight boundary conditions, the energy 
(Eshelby, 1956) that 

r 2 - 

[ £T° u re + Tre U 0e - <r re Ur - x r0e u° ] r=( b d6 = 0 
•*0 

must also be satisfied. The final result for the effective 
shear modulus is given by the quadratic equation 


(31) 


(32) 


, and rib. 


(33) 


(34) 


(35) 


(36) 


i are to be 
requirement 


(37) 


transverse 
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+ (^12/Mib) B + D - 0 

where 

A = 3V f (l-V f ) 2 [n f /|i Bl -l][fif/ji»+nf] + 

[ + VfVm ~ ( (Hf/Umlym ~ T Jf)Vf] x 

" ((Hr/nJVm +1)1. 


(38) 


(39) 


B = -6Vf ( 1-Vf ) + [(Mf/Mm^m + (Hf/Mm-l) V f + 

[(Tj m -l)(|Jf/Mm + Df) - 2Vf(() u c /nJv m - T) f )l + 

Vf ( T)m + 1 ) [ Mf /Mm- 1 ] ( Hf /Mm + Vf + ( (Hf/Mm^m ~ Vf )Vf ] , I 


D = 3V f (l-V f ) 2 [Mf/Mm-lHHf/Mm + nf] + [(Mf/Mo.)T) B + ( Mr /**m" 1 ) v f + H x 
[yLf/Hm + Vf + ((flf/MmJTJn. “ 


(41) 


The quadratic equation, Eq. 38, will have one negative root and one 
positive root, which is the effective transverse shear modulus, p i2 - 

The positive root of Eq. 38, plotted in Fig. 1, lies between both 
sets of Hashin bounds which demonstrates the validity of the solution. The 
one assumption made in the derivation of the Christensen-Lo model is that 
all fibers interact similarly with one another, and thus the outer phase 
can be modeled as an equivalent composite material. Even though this is a 
sound physical model, it is still an assumption that is required to solve 
the problem, and should be stated. Nevertheless, we shall refer to their 
result as the exact elasticity solution. The difficulty in using the 
three-phase model is that, like the Hashin model, it requires solving an 
elasticity problem. This makes it very difficult, if not impossible, to 
extend the method into the inelastic domain. In the self-consistent 
methods, the problem is formulated using the strains in tensor notation, 
which lends itself nicely to including inelastic strains. 


Self-Consistent Method for Non-Dilute Concentrations (NSC) 

In order to extend the basic self-consistent method to non-dilute 
concentrations of inclusions, Hill (1965) and Budiansky (1965) used the 
idea of an equivalent matrix surrounding the inclusion. This equivalent or 
homogeneous matrix was assumed to have the material properties of the 
overall composite, and thus accounts for the effect of the inclusions 
interacting with one another. 

The final form of Hill and Budiansky for predicting the effective 
material properties can be easily derived by starting with the results 
obtained in the DSC method. First, recall Eq. 4, and modify the D m tensor 
to be the D tensor which represents the effective material properties 
(remember, the fiber is now assumed to be embedded in a homogeneous 
composite material D, and not into the matrix D m ) . 
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D (e° + e c - c T ) = Df (e° + e c ) 


(42) 


Preceding in the same manner as with the DSC method, as outlined in Eqs. 5 
through 8, the final strain concentration tensor will be 

T = [I + S [D' 1 (D f - D)]]" 1 . (43) 

This can now be used in Eq. 13 to give 

D = D m + Vf (Df - DJ [I + S [D- 1 (Df - D)]]"\ (44) 

which is the standard _form of the NSC method. It should be noted also that 
the Eshelby tensor, S, is now a function of the homogenized composite 
phase. If the fibers are arranged In a_hexagonal array, the composite will 
be transversely isotropic and the S tensor will have the following 
components (Walker et al. 1989) 


S = 


5Dj i+Dj 2 3D 12 -D lt D 13 

8D| i QT>i \ 2 Dj i 

3Dj 2~Pi i 5Di i+Di 2 Pi 3 

8Dj i 8Dj i 2Dj i 

0 0 0 

0 0 0 

0 0 0 


0 0 

0 0 

0 0 

0.5 0 

0 0.5 


0 

0 

0 

0 

0 


0 


0 


0 


0 


3Dj 2 
4D lt 


(45) 


While Eq. 44 is not in the same form as the equations given by Hill 
(1965) and Budiansky (1965), it can be rearranged and manipulated to match 
their solutions. Nemat-Nasser (1981) has also derived an alternate form of 
Eq. 44, I.e. 

D = D m + Vf (Df - D m ) (D - Df)' 1 D [(D - Df)' 1 D - S ]"\ (46) 

where his notation has been changed to match the notation in this paper. 

Equation 46 can be manipulated as follows to produce Eq. 44. 

D = D m + Vf (D f - DJ (D - Df)' 1 [(D - Df)" 1 - S D" 1 ]' 1 , 

D = D m + Vf (Df - DJ [I - S D’ 1 (D - Df)]' 1 , 
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D = D m + V f (Df - D m ) [I + S [D-i (D f - D)]]' 1 . (Same as 44) 

Walker et al . (1989) have also developed an expression for the 
elastic effective material properties independent of the above methods. 
But as will be seen, their method is also another form of the basic NSC 
method. Their equation for the elastic effective material properties of a 
composite is given by Eq. 165 in Walker et al. (1989); it is 


D = V f D f [I + S [ D -1 (D f - D)J] 1 + 

+ ( 1— Vf) D m p + S [ D 1 (D m - D) ]] 


This form is difficult to manipulate to give Eq. 44. But if one returns to 
Eq. 114 in their report from which Eq. 47 was derived, an alternate 
derivation will give Eq. 44. Equation 114 of Walker - without the 
inelastic strain (eigenstrain) terms - is 


L&° = V f D f Ae T (f) + V m D m Ae T (m) , (48) 

which is the rule of mixtures for the stresses in an elastic composite. 
Furthermore, one observes that Eq. 48 is the incremental form of Eq. 9. 
Therefore, following the same method of derivation, Eq. 48 becomes Eq. 13, 

D = D m + V f (D f - DJ T, (13) 


where T is replaced by T. It Is interesting to note that Eq. 13 is also 
the same as Eq. 87 in Walker et al . (1989). The elastic portion of that 
equation is given by 


A<r° = D^Ac 0 - 6D(r) Ae T (r) dv, (49) 

where V is the total volume of the composite subvolume, and 5D is equal to 
0 in the matrix region and equal to Df-D m in the fiber region. Rewriting 
Eq. 49 in the standard notation used in this paper, and evaluating the 
integral gives 

A(T° = DmAc° - 0 - Vf [D m Ae T (f ) - D f Ae T (f)]. (50) 


This can be rearranged to give 


D Ae° = DmAe 0 + V f (D f - D m ) Ae T (f). 


(51) 
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This too is the same as Eq. 13, except that it is in an incremental form. 
Thus both Eqs. 87 and 114 of Walker’s are the same, and are also the same 
as Eq. 13 in this report. The only other term that is needed is the strain 
concentration tensor T. Again turning to the work of Walker et al. (1989), 
they derive this tensor for the inelastic case though extensive equation 
manipulation using a Green’s function formulation. The end result, given 
as Eq. 125 of their report, produces a T tensor identical to that in Eq. 
43 of this report. In summary, the elastic portion of their self- 
consistent method is the same as the standard NSC method. 

As expected, both Eq. 44 (the basic NSC equation) and Eq. 47 (the 
Walker solution) give the same results, but Eq. 44 converges much more 
rapidly. The results of the NSC method are shown in Figs. 1-6. In all 
cases, except for E 33 , the NSC results deviated from the elasticity 
solutions. This can be partially attributed to the assumption of constant 
strain In the Inclusions, which is not true at high concentrations of 
Inclusions. 

Useful insight can be acquired by comparing Budiansky’s (1965) 
equations with Eq. 44 for the bulk and shear moduli of a composite where 
the inclusions are spherical voids. Budiansky derived elastic solutions 
which show that if the bulk and shear moduli of the inclusions are zero 
valued (voids), and the matrix bulk modulus is infinite (incompressible), 
then the shear and bulk moduli of the composite will be given by 


V = 
K = 


r 3d-2Vf) L 

L 3-Vf J Mm ’ 

r 4( l-2V f ) ( 1-Vf ) 
[ V f (3-V f ) 


(52) 

(53) 


Equation 44, which Is in tensor form, can be rearranged to give Eq. 52, 
since there are no off-diagonal shear terms In the D matrix. However, Eq. 
53 cannot be obtained because the bulk modulus K is a combination of two 
or more D components, and can therefore only be solved numerically. But 
when Eq. 44 Is solved by successive substitution with the conditions = 
Kf = 0 and Km = oo, the solution diverges from the correct answer as given 
in Eqs. 52 and 53. One reason for this can be seen if D f = 0 Is substi- 
tuted into Eq. 44 to give 

D = D m - V f D m [I + S m [D-i (-D)]]"\ (54) 


or simplifying, 


D = D m [I - V f [I - Sj-»]. (55) 

The [I-S] term is singular whenever v = 0.5 is used to calculate the 
Eshelby tensor, S m . Also, D m is very large since K m = oo. These two 
conditions make it extremely hard to converge to the correct answer 
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numerically. Another numerical difficultly that arises in solving Eq. 54 
is when the Poisson ratio, v , is less then 0.2, which can easily happen 
during the convergence process. If v < 0.2 or v > 1.0, then the Eshelby 
tensor will have negative terms in the main diagonal of the matrix. This 
makes all direct substitution solution techniques diverge. 

The above two difficulties are not directly apparent for the case of 
fiber inclusions because of the complexity of the equations, but they may 
be one cause for the slow convergence. It should be noted that the D 
matrix does become symmetrical after the Iteration process has fully 
converged, thereby satisfying the requirement of a symmetric stiffness 
matrix. 


Revised Mori -Tanaka using the NSC Method 

Recall, in the derivation of the Mori and Tanaka (1973) method a new 
strain concentration matrix, G, was introduced which relates the strain in 
the inclusion to the strain in the matrix (e^ = G c m ) . This was then used 
to derive the following relation 

D = D ra + V f (D f - DJ A, (56) 
where A = G [V m I + VfG] 

At this point, Mori-Tanaka method sets G = T using ad hoc arguments that 
they advocate are a consequence of considering the non-dllute case. 
However, the actual form of T Is still derived from the dilute case. This 
can be shown by deriving the strain concentration matrix, G, In terms of 
the non-dllute condition. From Eq. 47 (or Eqs. 125 and 142 of Walker et 
al., 1989) the strain in the fiber and matrix for the non-dilute case can 


be written as 

e f = [i + S [D" 1 (Df - D)]]" 1 e° = T f c° (57) 

e m = [I + S [D- 1 (D m - D)]]" 1 c° = T m e° (58) 

Substituting Eqs. 57 and 58 into e f = G c m gives 

Tf e° = G T m c°, (59) 

and rearranging gives 

G = T f T‘\ (60) 


Looking at the limiting case as V f 0, then D -> D m and S -» S m , and G 
becomes 
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S T. 


(61) 


G -> G m = [I + S m [D^ 1 (Df - DJ3]' 1 

However, for the other limiting case as V f -) 1, then D D r and S S f , 
and G becomes 


G -> Gf = I + Sf [Df 1 (D ro - D f )] * T, (62) 

which is not the same T as implied in the original Mori-Tanaka method. 
Therefore, the Mori-Tanaka method cannot assume G = T for all cases of V r . 

It Is interesting to note that when G, as derived in Eq. 60, Is used 
in Eq. 56, the results for the effective material properties, D, are 
identical to the results from the NSC method. This should not be too 
surprising however, since the new G was derived using the same basic 
equations (Eqs. 47) for the matrix and fiber strain. The reduction of Eq. 
56 to the NSC equation is accomplished by substituting G back into Eq. 56 
giving 


A = Tf t; 1 [V m I + Vf Tf T; 1 ]" 1 (63) 

t 

which can be simplified to 

A = Tf [V m T m + V f Tf] -1 (64) 

The T terms in the brackets can be replaced by Eqs. 57 and 58, giving 

A = Tf e° [V ra e m + V f e f ] 1 (65) 


This reduces to 

A = Tf = [I + S [D-i (Df - D)]]" 1 (66) 

i 

because e° = V m e m + V f e f , which is the same as the NSC method. Therefore, 
this modified Mori-Tanaka method is just the NSC method. 


Summary and Conclusion 

Various self-consistent and elasticity solutions have been presented 
and compared with the method developed by Walker et al. (1989) for 
predicting the effective material properties. It was found that the Walker 
et al. method for the elastic material properties is identical to the 
self-consistent method for non-dilute concentrations (NSC) presented by 
Hill (1965) and Budiansky (1965), and later by Nemat-Nasser (1981). Also, 
the solutions for the NSC method are shown to deviate from the elastic 
solutions of Hashin (1965, 1972) and Christensen and Lo (1979), and the 
Mori-Tanaka (1973) self-consistent method, when the fiber volume ratio V f 
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exceeded about one-half. However, for practical metal matrix composites 
where the fiber volume ratio is generally less then one-half, all methods 
give adequate results. 

Because the NSC method does not stay within the accepted bounds for 
the effective elastic material properties, it should only be used with 
great care and caution. A second detriment to using the NSC method is the 
iterative nature of the solution. The convergence is slow using successive 
substitution, and would be difficult to implement into a finite element 
program where this method would be used at every element node point. It 
might be possible to use a higher order solution technique to improve on 
the convergence, but that was not investigated in this report. 

In order to stay within the elastic solution bounds, and to minimize 
the solution time, it is recommended that a non-iterative approach such as 
the Mori-Tanaka or Chrlstensen-Lo methods be examined for extension to 
include the inelastic material properties. The Christensen-Lo method would 
be the first choice because it is an exact elasticity solution, but it 
will be the more difficult of the two to extend into the inelastic region. 
The other choice, the Mori-Tanaka method, should be easier to extend 
because it is already in tensor form, but there are still unanswered 
questions about its theoretical validity. 
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Fig . 1 . Transverse Shear Modulus for a Glass/Polyester Fber Composite . 
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Fig . 2 . Longitudinal Shear Modulus for a Glass/Polyester Fiber Composite 
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Fig . 6 . Comparison of both Sets of Ha shin Bounds for a Glass/Polyester 
Fiber Composite with the non-cSute self-consistent method . 
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